OPTIMALITY OF GENERALIZED BERNSTEIN OPERATORS 



J. M. ALDAZ AND H. RENDER 

Abstract. We show that a certain optimahty property of the classical Bernstein op- 
erator also holds, when suitably reinterpreted, for generalized Bernstein operators on 
extended Chebyshev systems. 



1. Introduction 

It is well known that the classical Bernstein operator i?„ : C [0, 1] — t- Span{l, x, . . . ,x^}, 
defined by 



n—k 



k=0 

has good shape preserving properties, a fact that explains its usefulness in computer 
assisted geometric design. But the uniform convergence of the polynomials B^f to / 
can be slow. Thus, it is natural to enquire whether one can have operators from the 
continuous functions to the space of polynomials, with good shape preserving and fast 
approximation properties simultaneously. However, in a certain sense it is impossible to 
do better than by using Bernstein operators. Within some natural classes of polynomial 
valued operators, Bernstein operators approximate convex functions in an optimal way. 
In [5j, Berens and DeVore consider operators L„ : C [0, 1] — )■ Span{l,x, . . . fixing 1 

and X, and such that (Lnf)^^^ > if f^^'^ > for j = 0, 1, 2, ■ ■ ■ ,n. They show that for 
every n, Lnx"^ > Bnx"^, and if for some t G (0,1), Lnt"^ = Bnt'^, then L„ = Since 
BnX^ > x^, Ln performs strictly worse than Bn (unless L„ = Bn) at least for x^. The 
result of Berens and DeVore p. 214] is extended to a wider class of operators Ln in [6l 
Theorem 1]. While the result is stated there for x^, the argument works for an arbitrary 
convex function (using Jensen's inequality instead of Holder's inequality). The more 
general result, for convex, can be found in [71 Theorem 2]. A further extension appears 
in ^12i Theorem 2]. 

Recently, there has been an increasing interest concerning generalized Bernstein bases in 
the framework of extended Chebyshev spaces (in particular, spaces that contain functions 
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such as sinx and cosx), cf. for instance, [8], |9], ^0\, [H], [17], [H], [I9]. The existence 
of generahzed Bernstein bases raises the question of the possible existence of associated 
generahzed Bernstein operators. This topic is studied in [1], for exponential polynomials, 
and in [2], [3], [16], [20], for extended Chebyshev spaces. Here we show that the operators 
considered in [1], [2], [3], [16], and [20], share the preceding optimality property of the 
classical Bernstein operator, under a suitable generalized notion of convexity. 

Additionally, we give a selfcontained and streamlined presentation of the generalized 
Bernstein operators introduced in [2T] by S. Morigi and M. Neamtu. These operators are 
used here to exhibit a sequence of generalized Bernstein operators converging strongly to 
the identity, preserving (1, e^)- convexity, and failing to preserve standard convexity. 

2. Definitions and Optimality Results 

For each pair n, k of nonnegative integers, let fin,k be a positive measure. Denote by 
Pn,k{x) '■= (^)x'^ (1 — the elements of the Bernstein basis for the space of polyno- 

mials. Let Cn be the class of positive linear operators L„ : C [0, 1] — t- Span{l, x, . . . , x'^} 
fixing 1 and x (so preserves afiine functions) and defined by 

n „1 

(2) Lnf{x):=^\n,k{f)Pn,k{x), whcrc K,k{f) '■= I fdi^n,k- 

k=0 -^0 

Denote by Sx the Dirac delta probability, which assigns mass 1 to the singleton {x}, and 
observe that Bn G Cn is obtained by taking /i„ ^ = 6k/n in ([2]), so Xn,k{f) = f{k/n). Recall 
that if (f) is convex on [0, 1], then > on that interval. From [HI Theorem 1] and the 
generalization in [7, Theorem 2], it is known that L„0 > for every convex function 
(p and every L„ G £„, and furthermore, L„ = if Lnip{t) = Bnip{t) for some strictly 
convex function ijj and some t G (0, 1). In this sense the operator Bn approximates convex 
functions optimally within the class £„. We emphasize that this is not an asymptotic 
optimahty result, where n is required to approach oo. Instead, best approximation holds 
for each fixed n. 

Observe that by uniqueness of the coefficients of the basis functions, the condition 
L„l = 1 entails that each /i^^ is a probability measure, while L„x = x entails that the 
expectation of x with respect to fin,k is Jq xdfin,kix) = k/n. If the assumption L„x = x is 
dropped, then the resulting operator may fail to preserve convexity, and may approximate 
some convex functions better than Bn (while in other cases the approximation will be 
worse, notably for x itself). 

Example 1. Let 1 < j < n. In [2", Proposition 11] and in [201 Example 5.3], a generalized 
Bernstein operator Bn,j fixing 1 and x^ , is defined by 

n 

BnJ{x) = J2f 
k=0 



'k{k-l)---{k-j + 1) 
n{n — 1) • • • (n — j + 1) 




X il - X 



\n—k 
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In particular, if j = 1 this is the classical Bernstein operator. However, for j > 1 and 
all k = 1,2,..., BnjX^ < BnX^ on (0,1), since the operator samples the strictly 
increasing functions at nodes lower than those of (except the first and the last 
nodes, which are and 1 in both cases). It is shown in [21 Example 17] that if 1 < A; < j, 
then Bnjx'' < x'^, while if j < /c, then Bn,jX^ > x^ . Hence, on (0, 1), BnX^ > Bn,jX^ > x^ 
when i < k. Of course, Bnj does not preserve convexity, since 

i3„j(l -x)= Bnj{l) - Bnj{x) > 1 - B.^{x) = l-x 

whenever x G (0, 1). But Bn,j does preserve (1, a:-') -convexity, defined below. 

The way to recover the results from [6] and [7] is thus to deal with a generalized notion 
of convexity. We give next the relevant definitions. An extended Chebyshev space Un of 
dimension n + 1 over the interval [a, 6] is an n + 1 dimensional subspace of C"'{[a,b]) 
such that each function f & Un not vanishing identically, has at most n zeros in [a, b] , 
counting multiplicities. It is well-known that extended Chebyshev spaces possess non- 
negative Bernstein bases, i.e. collections of non-negative functions Pn,k,k = 0,...,n, in 
Un, such that each pn,k has a zero of order at a and a zero of order n — k ai b, for 
k = 0,...,n. Now let us select two functions /o, /i G Un, such that /o > and /i//o is 
strictly increasing (these functions play the role of 1 and x in the classical case). It is 
sometimes possible to define a generalized Bernstein operator Bn '■ C [a, b] — )■ Un fixing /q 
and /i, by first suitably choosing nodes tn,o, tn,n ^ [^5 b] such that tnfi = a and tn,n = b, 
second, by selecting appropriate weights anfi, > 0, and finally, by setting 

n 

(3) Bn {f)=J2f 

n,k ) (-^n,kPn,k- 

k=0 

Criteria for the existence of such Bn can be found in |2] , [3] , |20] , and in the specific case 
of exponential polynomials, in pQ and (TB]. In this section, we shall assume that Bn exists 
(fixing the given functions), and will show that it has an extremal property analogous to 
that of Bn- For simplicity in the notation, we relabel the basis functions so that they 
already incorporate the weights an,k- Thus, becomes 

n 

(4) Bn{f) = J2f(tn,k)Pn,k. 

k=0 

The functions /o and fi are always assumed to satisfy the conditions "/o > 0" and "/i//o 
is strictly increasing" . 

Remark 2. Given an extended Chebyshev space [/„, it is always possible to select some 
pair /o and /i in Un for which a generalized Bernstein operator can be defined (cf. [21 
Corollary 8]). However, for some choices of /o and /i, Bn may fail to exist. For instance, 
it is well known that if < 6 < 2n, then U3 := Span{l, x, cosx, sinx} on the interval [0, b] 
is an extended Chebyshev space. But there is no generalized Bernstein operator Bn fixing 
1 and a; if 6 = 5 (for instance) while such Bn does exist if 6 = 4, cf. [21 Theorem 25]. 
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Remark 3. The condition on the zeros of pn,k entails that, hke the classical Bernstein 
operator, the operators given by interpolate functions at the endpoints of the interval 
[a, b]: Since /q is strictly positive, from fo{a) = i3„ (/o) (a) = /o (a) an,oPn,o('2) "we conclude 
that anoPno{ci) = 1- Likewise, OinnVnnip) = 1, so we always have i3„ (/) (a) = /(a) and 
B^{f)(b) = f{b). 

Convexity of can be defined by saying that if we interpolate between (j){x) and (j){y) 
using an affine function h, then on [x, y] the graph of </> lies below the graph of h. If in- 
stead of affine functions we use (/o, /i)-affine functions, that is, functions in Span{/o, /i}, 
we obtain the corresponding notion of (/o, fi)- convexity. Thus, ordinary convexity corre- 
sponds to (1, x)-convexity. According to [H], p. 376, this generalized notion of convexity 
was introduced in 1926 by Hopf, and was later extensively developed by Popoviciu, spe- 
cially in the context of Chebyshev spaces. Strict (/o, /i)-convexity is defined analogously 
to strict convexity. It is shown in [21 Theorem 22] (cf. also ^20^ Proposition 4.14]) that a 
generalized Bernstein operator i3„ fixing /q and /i preserves (/o, /i)- convexity. 

We would like to stress the point that i3„ will in general fail to preserve standard 
convexity (recall example [T]). It will be shown in Section 3 that likewise, the standard 
Bernstein operator does not in general preserve strengthened forms of convexity, such 
as (1, e^)-convexity for increasing functions. 

The following characterization of (/o, /i)-convexity, due to M. Bessenyei and Z. Pales 
(cf. [H Theorem 5]), helps to understand its meaning. 

Theorem 4. Let I := (/i//o)([Oj &])■ Then G C[a,6] is {Jq, fi)-convex if and only if 
(0//o) ° (/i//o)^^ £ C (/) is convex in the standard sense. 

In particular, if /o = 1, the relationship between convexity and (1, /i)-convexity is given 
by a simple change of variables (determined by /i). For instance, returning to Example [H 
and to the functions considered there, we see that ii I < k < j, then x'' is (1, x-')-concave, 
while ii j < k < n, then x^ is (1, x-^)-convex. 

Before we prove the announced optimality property of generalized Bernstein operators, 
we need the following lemma. It is entirely analogous to Theorem m save that strict 
convexity replaces convexity throughout. The proof is also essentially the same, and thus 
we omit it. Basically, all one needs to do is to use strict inequalities > instead of >, at 
the appropriate places. 

Lemma 5. Let I := (/i//o)([fl, &])■ Thenip G C[a,6] is strictly {f fi)- convex if and only 
if ii'/fo) o {fi/foY^ e C (/) is strictly convex. 

Theorem 6. Let Un C C" [a, b] he an extended Chebyshev space containing the functions 
fo and fi, with fo > and /i//o strictly increasing. Suppose there exists a generalized 
Bernstein operator Bn '■ C [a, b] — > Un fixing /o and fi . Let the operator Ln : C [a, b] ^ Un, 
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subject to Lnfo = fo and Lnfi = fi, be defined by 

n 

(5) Ln if) = ^ Xn,kif)Pn,k, 

k=0 

where Pn,k has the same meaning as in and Xn,k{f) is obtained from f and the positive 
measure fin,k via 

(6) Xn,k{f) ■= foitn,k) / fdHn,k- 

J a 

Then, for every {fo, fi) -convex function (j), we have < < Lncj). Moreover, if for 
some t G (a, 6) and some strictly {fo, fx) -convex function ip we have Bni'it) = Lni'it), 
then Bn = Ln- 

Proof. First we make explicit what the assumptions in the theorem entail about the 
measures iJin,k- Since 

n n 

Lnfo = ^ Xn,k{fo)Pn,k = fo = ^n/o = ^ /o (^n,fc)Pn,fc, 
fc=0 k=0 

by uniqueness of the coefficients of the basis functions, we must have fo {tn,k) = Xn,k{fo) 
for all pairs n,k. But then 1 = fodfin,k by IQ, so admissible measures fin,k niust be 
such that dPn,k '■= fodfJ,n,k defines a probability Pn,k- With this notation, the second 
condition = fi leads to 



Xn,k{fl) — fo (tn,k 



I ^dPn,k, 
J a JO 



SO again by uniqueness of the coefficients of the basis functions, we have 

fo {tn,k) [ ^dPn,k = A„,fc(/l) = fl (tn.k) = fo (tn.k) ^1 / • 

J a Jo Jo [tn,k) 

Thus, 

(7) TTr\= fr'"'-'^ 

Jo [l^n,k) J a Jo 

for all pairs n, k. 

Next, suppose is (/o, /i)-convex on [a, h]. Given h : X Y and a measure v on X, 
the pushforward to Y of z/, using h, is denoted by h^v] recall that h^p^A) := z/(/i^^(A)). 
Applying the change of variables formula twice (cf., for instance, [TBI Theorem C, p. 163]), 
Theorem m above, Jensen's inequality, and (I7j), we obtain 



K,k{<P) = fo{tn,k) I 4>{x)dlJn,k{^) = fo {tn,k) —{x)dPn,k{x) 

Ja JO 
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(9) = /O {tn,k) I ^ O ( ^1 i^)d (fy ] Pn,k{^) 



/o' 



j^{la,b]) \foJ \fo 



(10) >/o(tn..)^o( A) I / :,d{^] P„,fc(x) 

Jo V Jo , 



Jo {tn,k) -r° [-r 
Jo V/o 




(11) = /o {tn,k) ^°{-rJ (I i-r] ix)dPn,k{x 

(12) = /O {tn,k) I O j (t^^^"'^) j = '^^^"''^-^ 

Since the basis functions are non-negative, we obtain i3„0 < L„0 just by adding up. And 
the fact that < i3„0 for all (/o, /i)-convex functions is proven in ^ Theorem 15]. 

Assume next that for some t G (a, b) and some strictly (/o, /i)-convex function ip we 
have Bnipit) = Lnipit). We must show that Pn,k = St„ ^ ^"^^ QSi.ch. pair tt,, k. Set = in 
(11- (fH. Since 



A„,fc(V')Pn,fc(^) = ^nV'(^) = 'B„V(^) = (^n,fe)Pn,fc(t), 
fc=0 fc=0 

and for all n, /c, we have Xn,k{'4') ^ '0 (^n,fc) (by fl8|)- f|T2|) ) from pn,k > on (a, b) we conclude 
that Xn,k{'^) = "0 (^n,fc) for all n, k. Thus, with (p = ip we have equality in ffTOj) . By Lemma 

[5l the function j^o is strictly convex. So by the equality case in Jensen's inequality 

we conclude that the function x is constant a.e. with respect to Pn,k- This entails 

that the latter measure is a Dirac delta, and since /i//o is strictly increasing, P„ fc must 
also be a Dirac delta. It now follows from ([7]) that Pn^k = St„ k ■ ^ 

3. On the generalized Bernstein Operator of Morigi and Neamtu 

There are known examples of chains of extended Chebyshev spaces, for which their 
associated generalized Bernstein operators exist and converge strongly to the identity. 
An instance is given by the operators Bnj described in Example [T] (cf. [21 Proposition 
11] or |2ni Theorem 6.1] for the convergence assertion). As was observed in Example [H 
the operators Bnj preserve (1, a;-') -convexity, and fail to preserve ordinary convexity when 
J > 1- 

Nevertheless, these operators are somewhat degenerate in that the sequences of nodes 
are not strictly increasing, due to the fact that for j > 1, /i(x) = has a zero of order at 
least two at zero, and in particular, f[{0) = 0. Thus, f[ is not strictly positive on [0, 1]. 
We shall use the generalized Bernstein operators introduced by Morigi and Neamtu in 
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[2T| to present a sequence of operators that i) converges to the identity, ii) the operators 
preserve (1, e^)-convexity, iii) they do not preserve ordinary convexity, iv) they are defined 
via strictly increasing sequences of nodes (in fact, the nodes are equidistributed) , and v) 
they fix /o = 1 and /i(x)//o(x) = /i(x) = e^. Obviously we have /( > everywhere. 

Our presentation here of the Morigi and Neamtu operators is selfcontained and simplifies 
some of the arguments from [21]. At the beginning we will consider a setting more general 
than we really need (from the shape preservation perspective) allowing complex values. 
We start by defining these operators. 

Given two different complex numbers /io,Aii, set u;„ := ^ (/ii — /io) 7^ 0, and for j = 
0, n, let 

(13) Xj:=fio+jUn. 
Thus, Ao = yUo and A„ = /ii. Define 

(14) ^„ (x) := ^ _ and (x) := 1 - ^„ (x) = ^ _ • 

Then (fn)'' is a linear combination of the exponential functions e'^"^^, where j = 0, k, 
and {ipn)"'~^ is a linear combination of the functions e'^"-'^, where j = 0, n — k,. There- 
fore, 

(15) p.,, (x) := e-^"t(^-) . e^o(--) Q (^„ {x)f (^„ (x))"-'= 

is a linear combination of the exponential functions e^^^', j = 0, ...,n. It follows that p^.fc 
is contained in the vector space 

(16) c„:={/,C(«.C):(^-A,)...(A_A„)/ = 

Since has a simple zero at x = a and ipn a simple zero at b, pn,k has a zero order k at 
a and a zero of order n — k aih. 

It is well known that if Aq, A„ are real, then := {f G Un '■ f is real-valued} is an 
extended Chebyshev space over any bounded closed interval [a, b] (we mention that when 
deahng with complex values Aq, A„, restrictions must be imposed on the length of [a, b] 
in order to ensure that f/„ is an extended Chebyshev space, cf. for instance [1), Theorem 
19]). Furthermore, {pn,k : = 0, . . . , n} is a non-negative Bernstein basis. 

Next we specify nodes and weights. Let t^^k '■= o. + ^ {b — a) for = 0, n, and let 

n 

(17) 13n{f){x):=J2fitn,k)Pn,k{x). 

k=0 

Proposition 7. The operator Bn fixes e'^"^ and e^^^ . 
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Proof. Recall from ( fT3l) that /io = Aq and /xi = A„. Using ( IT71) . ( |T5l) . and ( Tl4l) . we obtain 

i3n (e'°^) = e^»^ r^") (1 - {x)r-' = e^°^ (x) + 1 - (x))" = e^°^ 

^— n \ / 



k=0 

Furthermore, 



fc=0 ^ ^ 



,A„agAo(x-a) (e-n(fc-«)y,^ (x) + 1 - y5„ (x))" 

But now (IT^ implies that 

^g.;„(6-a) _ 1) (x) + 1 = e""(^-") (1 - e"("-'')) (x) + 1 

SO from the preceding equalities we obtain 



□ 



For simplicity, next we take [a, b] to be [0, 1], /io = 0, and fii = 1. Then the Morigi and 
Neamtu operator reduces to 



k=0 V / \ / 



18) ^nfix):=yj (x))'= (1 - ^„ (x)) 



n—k 



where by specialization of f fT4|) we have 

Observe that since e^^" is convex and increasing, so is <y9„. Furthermore, <^n (0) = and 
ipn (1) = 1- Recalling that Bn denotes the classical Bernstein operator, given by ([1]), we 
see that 

(19) BJ{x) = BJ{^^{x)). 

From the preceding expression it is immediate that Bn does not in general preserve the 
convexity of a (decreasing) convex function. Let h{x) = — x. Since Bnh = h, we have 

Bnh (x) = Bnh {ipn (x)) = -ipn (x) . 
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But (fn is strictly convex, so the image under Bn of the affine (hence convex) function h 
is the strictly concave function — v^n- This shows that the following assertion, contained 
in [201 Remark 4.15 (1)], is not correct: If /o = 1, /i is convex, and > everywhere on 
[a, 6], then the image of a decreasing convex function under i3„ is convex (this is claimed 
to follow from [20i Proposition 4.14], but it does not; we note that Proposition 4.14 from 
[20] is correct, cf. [21 Theorem 22] together with Theorem 4 above for a stronger result). 

It is shown in [211 Theorem 3.3] that the operators Bn converge strongly to the identity 
(more general or related results can be found in [U Theorem 23] and [SUl Theorem 6.1]). 
However, for the special case we are considering, it is easy to give a direct argument, so 
in order to make this section as selfcontained as possible, we do this next. 

Let g{x) := x. Since Bn fixes 1 and e^, by Korovkin's Theorem it suffices to prove that 
Bng converges uniformly to g on [0, 1]. By f|T9|) . Bng = ^Pn, so it is enough to show that 
ipn converges uniformly to g. Since (pn is convex and increasing, its derivative achieves 
its maximum value over [0, 1] at 1. Now let gn be the affine function with slope v'^(l) 
passing thorough the point (1, 1). By convexity, gn 'Pn g, so to obtain the uniform 
convergence it suffices to show that lim„_^oo V^n(l) = 1- But this follows immediately by 
writing explicitly 93^(1), and then using L'Hospital rule, or the Taylor expansion for e^^". 

We have seen that Bn does not preserve convexity. It is easy to check that Bn does not 
preserve (1, e^)-convexity either. Actually, both of these facts can be derived from the 
slightly more general result given next. 

Proposition 8. Let fo = l and fi be functions in an Extended Chebyshev Space E over 
the interval [a, h], where fi is assumed to be strictly increasing and strictly convex. Suppose 
there exists a generalized Bernstein operator Bn '■ C[a,b] — )■ E fixing fo and fi. Then the 
classical Bernstein operator Bn does not preserve {1, fi)- convexity, and the generalized 
operator Bn does not preserve convexity. 

Proof. For the first assertion, observe that while /i is trivially (1, /i)-convex, Bnfi is not. 
To see why, note that since Bnfi{a) = /i(a) and Bnfiip) = fi{b), the unique function 
in Span{l,/i} interpolating -B„/i at a and b is /i itself. If -B„/i were (1, /i)-convex, 
we would have -B„/i < /i on [a,b]. But /i is strictly convex, so by Jensen's inequality 
(including the equality condition) we have Bnfi{x) > fi{x) for all x G (a, 6). 

An analogous argument yields the second assertion. Set h{x) = x, so both h and —h 
are convex. Towards a contradiction, suppose that Bn preserves convexity. Then Bnh < h 
and Bn{—h) < —h, so Bnh = h. Fix x G {a,b), and note that / Bnf{x) is a positive 
linear functional on C[a,b], defined by a probability measure (since Bnl{x) = 1). Thus 
we can once more apply Jensen's inequality, in this case to the strictly concave function 
/f^, to conclude that 

X = h{x) = Bnh{x) = Bnif,-' O /i)(x) > /^^ (Sn/l(x)) = f^' (/l(x)) = X. 
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Alternatively, we can reach the same conclusion without using Jensen's inequality, by 
observing that if in addition to fixing 1 and /i, we have Bnh = h, then Bn fixes three lin- 
early independent functions, and thus it must be the identity operator (by [20, Proposition 
3.7]). So again we have a contradiction. □ 

We finish with an easy Proposition, which sheds some light on the relationship between 
the different notions of convexity considered above. 

Proposition 9. Let fi be increasing and strictly convex on [a, b]. For increasing functions 
on [a,b], the condition of {1, fi)- convexity is strictly stronger than convexity. Over the 
decreasing functions on [a,b], convexity is strictly stronger than {1, fi)- convexity. 

Proof. Let / be increasing and (1, /i)-convex on [a, b]. Pick x,y E [a, b] with x < y. Since 
f{y) > f{^): the unique interpolant ip of f at x and y such that ip G Span{l,/i} must 
be of the form = cq + ci/i, with ci > 0. Thus ip is convex, so the unique interpolant 
(j) G Span{l,t} of ip (or /) at x and y satisfies (f){t) > i/j(t) > f(t) for all t G [x,y], 
whence / is convex. Also, h{t) := t is increasing and convex but not (1, /i)-convex, since 
h{fi^) = f{^ is strictly concave, cf. Lemma [51 

The assertion about decreasing functions is proved by the same type of argument, so 
we omit it. □ 



References 

[1] J. M. Aldaz, O. Kounchev, H. Render, Bernstein operators for exponential polynomials. Constr. 

Approx. 29 (2009), no. 3, 345-367. 
[2] J. M. Aldaz, O. Kounchev, H. Render, Shape preserving properties of Generalized Bernstein operators 

for extended Chebyshev systems. Numer. Math. 114 (2009) 1-25, DOI 10.1007/s00211-009-0248-0. 

Also available at the Math. ArXiv, since 12 May 2008. 
[3] J. M. Aldaz, O. Kounchev, H. Render, Bernstein operators for extended Chebyshev systems, submit- 
ted. Available at the Math. ArXiv since, 12 May 2008. 
[4] M. Bessenyei, Z. Pales, Hadamard-type inequalities for generalized convex functions, Math. Ineq. 

Appl. vol. 6, n. 3 (2003), 379-392. 
[5] H. Berens, R. DeVore, A characterization of Bernstein polynomials. Approximation theory. III (Proc. 

Conf., Univ. Texas, Austin, Tex., 1980), pp. 213-219, Academic Press, New York-London, 1980. 
[6] J. Bustamante, J. M. Quesada, On an extremal relation of Bernstein operators. J. Approx. Theory 

141 (2006), no. 2, 214-215. 
[7] Jesus de la Cal, Javier Carcamo, An extremal property of Bernstein operators. J. Approx. Theory 

146 (2007), no. 1, 87-90. 

[8] J. M. Carnicer, E. Mainar, J.M. Peiia, Critical Length for Design Purposes and Extended Chebyshev 

Spaces, Constr. Approx. 20 (2004), 55-71. 
[9] J.M. Carnicer, E. Mainar, J.M. Peiia, Shape preservation regions for six- dimensional space. Advances 
in Computational Mathematics 26 (2007), 121-136. 
[10] P. Costantini, Curve and surface construction using variable degree polynomial splines, Comput. 

Aided Geom. Design 17 (2000), 419-446, 
[11] P. Costantini, T. Lyche, C. Manni, On a class of weak Tchebycheff systems, Numer. Math. 101 
(2005), 333-354. 



OPTIMALITY OF GENERALIZED BERNSTEIN OPERATORS 



11 



[12] loan Gavrca, Mircca Ivan, An extremal property for a class of positive linear operators, J. of Approx. 

Theory, Vol. 162, 1 (2010) 6 9. 
[13] Paul R. Halmos, Measure Theory. D. Van Nostrand Company, Inc., New York, N. Y., 1950. xi+304 

pp. 

[14] S. Karlin, W. J. Studden, Tchebycheff Systems: with applications in Analysis and Statistics. John 

Wiley and Sons, 1966. 
[15] S. Karlin, Total positivity, Vol. 1, Stanford Univ. Press, Standford 1968. 

[16] O. Kounchev, H. Render, On the Bernstein operator of S. Morigi and M. Neamtu, Result. Math. 53, 
311-322. 

[17] E. Mainar, J.M. Pena, J. Sanchez-Reyes, Shape preserving alternatives to the rational Bezier model, 

Comput. Aided Geom. Design 14 (1997), 5-11, 
[18] M. Mazure, Bernstein bases in Muntz spaces, Numer. Algor. 22 (1999), 285 304. 
[19] M. Mazure, Chebyshev Spaces and Bernstein bases, Constr. Approx. 22 (2005), 347-363. 
[20] M. Mazure, Bernstein-type operators in Chebyshev Spaces, Nunicr. Algor. 52 (2009), 93-128. 
[21] S. Morigi, M. Neamtu, Some results for a class of generalized polynomials. Adv. Comput. Math. 12 

(2000), 133-149. 

J. M. Aldaz: Departamento de Matematicas, Universidad Autonoma de Madrid, Can- 
TOBLANCO 28049, Madrid, Spain. 

E-mail address: jesus.muiiarriz@uaiii.es 

H. Render: School of Mathematical Sciences, University College Dublin, Dublin 4, 
Ireland. 

E-mail address: renderOgmx . de 



